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Abstract 

By introducing a new entangled state representation, we show that the Laguerre-Gaussian 
(LG) mode is just the wave function of the common eigenvector of the orbital angular momentum 
and the total photon number operators of 2-d oscillator, which can be generated by 50:50 beam 
splitter with the phase difference <f) — n/2 between the reflected and transmitted fields. Based on 
this and using the Weyl ordering invariance under similar transforms, the Wigner representation 
of LG is directly obtained, which can be considered as the generalized Wigner transform of 
Hermite Gaussian modes. 

PACS: 03.65.-w-Quantum mechanics 

PACS: 42.50.-p-Quantum optics 

1 Introduction 

It has been known that a Laguerre-Gaussian (LG) beam of paraxial light has a well-defined orbital 
angular momentum [IH5], which is useful in studying quantum entanglement [6]. In Ref. [3] Nien- 
huis and Allen employed operator algebra to describe the Laguerre-Gaussian beam, and noticed 
that Laguerre-Gaussian modes are laser mode analog of the angular momentum eigenstates of the 
isotropic 2-d harmonic oscillator. In Ref. [7 Simon and Agarwal presented a phase-space description 
(the Wigner function) of the LG mode by exploiting the underlying phase-space symmetry. In this 
Letter we shall go a step further to show that LG mode is just the wave function of the common 
eigenvector |n, 1} of the orbital angular momentum operator and the total photon number operator 
of 2-d oscillator in the entangled state representation (ESR) . The ESR was constructed [HIE] based 
on the Einstein- Podolsky- Rosen quantum entanglement [TU]. It is shown that \n, I) can be generated 
by 50:50 beam splitter with the phase difference 4> = ir/2 between the reflected and transmitted 
fields. Then we use the Weyl ordering form of the Wigner operator and the Weyl ordering's co- 
variance under similar transformations to directly derive the Wigner representation of LG beams, 
which seems economical. The marginal distributions of Wigner function (WF) are also obtained 
by the entangled state representation. It is found that the amplitude of marginal distribution is 
just the eigenfunction of the fractional Fourier transform (FrFT). In addition, LG mode can also be 
considered as the generalized Wigner transform of Hermite Gaussian modes by using the Schmidt 
decomposition of the ESR. 



2 Eigenvector corresponding to Laguerre- Guassian mode 

In Ref. [3] the Bosonic operator algebra of the quantum harmonic oscillator is applied to the de- 
scription of Gaussian modes of a laser beam, i.e., a paraxial beam of light is described by operators' 
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eigenvector equations 



N \n, I) — n \n, I) , N = (a\ai + a\o,2j j 
L|n,/) = l|n,l), i = X 1 P 2 -X 2 P 1 , 



(1) 

since [N, L] — 0, where a\ and <Zj (i = 1,2) are Bose creation operator and annihilation operator; 
L and N are the orbital angular momentum operator and the total photon number operator of a 

paraxial beam of light, respectively. Using Xi = (a,i + ajj /y/2 and p = (a,i — a] J /(i\/2), and 



5<j, then 



L = i(a 2 ai ~ ffliffls)- (2) 
Here we search for the common eigenvector of (JV, L) in the entangled state representation. By 
introducing 

i i 

(3) 



which obey the commutative relation 

[Ah4] = 1, [A_,4 = l, 

[A + ,4=0, [A_,4=o, 

one can see 



TV = AT A 4 



A f _A_, L = A\A+ - aIa_. 



Now we introduce the entangled state representation in Fock space, 

\rf) = cxp |-i |r/| 2 + r]A\ - rf Al + A\Al j |00) , 



(4) 



(5) 



(6) 



here r\ = \r]\ e ltp = r\\ + ir)2, 1 00) is annihilated by A + and A_. It is not difficult to see that \rf) is the 
common eigenvector of operators (Xi — X 2 — Pi + P 2 , Pi + P 2 — -Xi — X 2 ) 

{X 1 -X 2 -Pi + P 2 )\r l }=2r h \r l }, 

(P + P 2 - X x - X 2 ) \rj) = 2 m \q) . (7) 
Using the normal ordering form of vacuum projector 

1 00) (00| = : cxp (-A\A+ - AlAJ) : , (8) 

(where : : denotes normal ordering) and the technique of integral within an ordered product (IWOP) 
of operators [HUE] we can prove the completeness relation and the orthonormal property of , 

d 2 r] 



Thus 1 77) is qualified to make up a new representation. It follows from ^ and ^ that 

a + \ V ) = ( V + Al) \ V ) , 4 \ V ) = (J- + ^j \ v ) , 



(9) 

(10) 
(11) 



which lead to (denote r — \rj\ for simplicity) 

1 



(AlA + + AlA-)\v)= [-r 2 -l-2 



drjdri* 
d 2 



r z 1 
2 2 V dr 2 ' r dr 



\n) 



Id Id 2 



dip 7 



\n) 



(AU + - AlAJ) \r,) = [v— - rf — \v) = \v) 



drj dr]* 



dip 



(12) 
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Projecting Eqs.([T]) onto the (rj\ representation and using ([3]), (ITU|) -(' |20j). one can obtain the following 
equations 

l(r,\n,l)=i-^-(r 1 \n,l), (13) 



and 



n (r)\ n, I) 



r 



1 (cP_ ld_ 1 d 2 

2 \ dr 2 r dr r2 dip 2 



(v\ n,l) 



(14) 



Eq. lflU)) indicates that (rj\ n, I) oc e . From the uniqueness of wave function, e ''^1^=0 = e l ' ¥ '| t ^ = 27r, 
we know Z = 0, ±1, ±2 • • • .So letting (rj\ n, I) — R(r)e~' dip and substituting it into (TUT) yields 



d 2 i? 1 dR 



/ 2 



(15) 



Introducing £ = r 2 such that 



Eq. (TTSJ) becomes 



— - U— — - 0— A?— 
~ Ve d£ ' ~d^ ~ It + ^d?' 



d 2 i? 1 dR 
Then make the variable transform in (fTT|) 



1 71 + 1 Z 2 

4 + ~~2£ 4^ 



i? = 0. 



A(0 = e-^w/vo 



one can obtain the equation for it(£), 



(16) 



(17) 



(18) 



du n - 



'-u = 0. 



(19) 



Ea. (|19[) is just a confluent hypergeometric equation whose solution is the associate Laguerre poly- 
nomials, £„!,(£), where n p = " 2 ^ , (n p = 0, 1, 2, • • • ) [T3]. Thus the wave function of \n, I) in (rj\ 
representation is given by 



(77| n,/)=C ie -^e-^ 2 rl'l^l(r 2 ), 



(20) 

where C\ is an integral constant. The right-hand side of Ea. (j20"|) is just the LG mode, so we reach 
the conclusion that the wave function of \n,l) in the entangled state representation is just the LG 
mode, i.e., the LG mode gets its new physical meaning in quantum optics. 

Next, we further derive the explicit expression of \n, I) . Using the completeness relation of (t/| (|9]) 
and (f2"0")l . we have 



\n,l) 



d 2 T) 



1 77) (77KO 



d 2 V -i r 2 I \ -Hip 121 t \l\ 1 2\ 

e 2 \r\)e ^r 1 l L% (r ). 



(21) 



Then noticing the relation between two- variable Hermite polynomial [14U15] and Laguerre polyno- 
mial, 

H m . n (77,77*) = ml (-irr,™-™^-™ (7777*) , (22) 



where m < n, and the generating function of H m ^ n (77, 77*) is 

Qm-\-n 



H rn ,n (x,y) 



dt m dt' n 



exp \-tt' + tx + t'y] 



(23) 
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as well as using the integral formula [TB] 

J ^ exp (C \z\ 2 +Hz + V z*) = ~e~f , Re (£) < 0, (24) 
we can reform Ea. (|2"Tj) as (without loss of the generality, setting I > and m p = [n + \l\]/2) 



(-lf'd f d 2 rj 



n p \ 



'-H np , m Av,V*)e- m \V) 



(-l)" p Ci d n » +m " 
np~\ dt n pdt' m p 



exp 



d 2 ri 
x / exp 

7T 



\r)f + (Al + t)r)+ (t'-Al)r) 
(-!)"" d d n " +m p 



t=t'=0 



|00) 



n p l dt n pdt' m p 



exp 



AW - tAl 



t=t'=0 



|00) 



(25) 



3 Generation of |n, /) by Beam Splitter 

Note Eq.© and 



Al = e^^e"^' 7 *, A[_ = e l ^ J *a\e 



J x = \ (a\a 2 + ajai) , |00) = e^ J * |00) 
thus Eq. (j2"5|) can be further put into the following form 

\n,l) = ^e^(a\) mp (a^ |00) 

\m p ,n p ) . 



(26) 



C\ n 



(27) 



It is easy to see that the normalized constant can be chosen as C\ — yfn p \/m p \, which further leads 
to 

\n,l) = e l % J * \m p ,n p ), (28) 

where J x can be expressed by angular momentum operators J+ = a\a 2 and J_ = a\a} 2 , J x = 
\ (J+ + J_) . J+, J- and J z — i ^ajai — a^^j make up a close SU(2) Lie algebra. 

On the other hand, the beam splitter is one of the few experimentally accessible devices that may 
act as an entangler. In fact, the role of a beam splitter operator [T7HT5] is expressed by 



B | 



exp 



(29) 



with the amplitude reflection and transmission coefficients T — cos|, R = sin|. The beam split- 
ter gives the phase difference <fi between the reflected and transmitted fields. Comparing Eg . (j2"9"]l 
with e l ^ J " leads us to choose 9 = ir/2 (corresponding to 50:50 beam splitter) and <\> = tt/2, thus 
B (tt/2, tt/2) is just equivalent to e l ^ Jx in form. This indicates that |n, 1} can be generated by acting a 
symmetric beam splitter with <p = 7r/2 on two independent input Fock states \m p , n p ) 



l Pl\ \ lb Pl2 ■ 



In addition, note that n = m p + n p , i.e., when the total number of input photons is n, so the output 
state becomes an (n + l)-dimensional entangled state [19| . 
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4 The Wigner representation 

As is well-known, the Wigner quasidistribution provides with a definite phase space distribution of 
quantum states and is very useful in quantum statistics and quantum optics. In this section, we 
evaluate the Wigner representation of |n, I) . According Ref. [3D], the Wigner representation of \n, I) 
is given by 

W\ n> i) = (n,l\Ax (x 1 ,p 1 )A 2 (x 2 ,p 2 ) \n,l) 
= (m p , n p \ e~ l ^ J *A 1 (xi,pi) 

x A 2 {x 2 ,p 2 )eJ% J * \m p ,n p ) , (30) 
where Ai (xi,pi) is the single- mode Wigner operator, whose Weyl ordering form [21] is 

A 1 (x 1 ,p 1 )= \d(p 1 -P 1 )S(x 1 -X 1 ) ;, (31) 

where the symbol; ; denotes Weyl ordering [32]. Note that the order of Bose operators a and a' 
within ; ; can be permitted. That is to say, even though [a, a'*] = 1, we can have [air \ — \a/a\. 
According to the covariance property of Weyl ordering under similar transformations |21j and 

v 2 

e^ J -X 2 e^ J ' = -^=(X 2 -P 1 ), 

e -^ J *P 2 e^ J * = -^(P 2 + Xt) , (32) 
V2 



we have 



e-^ J *A 1 (x 1 , Pl )A 2 (x 2 ,p 2 )e^ J * 

:,( Pi+X 2 \ x f X x -P 2 
= 6 \pi -= — )0 \ Xi 



V2 J V V2 

<l P z+ X i\ x { X 2 -P l 
xd [p 2 -= — ) S \x 2 



V2 J V V2 

. X X +p 2 pi - X 2 \ f X 2 +Pl P2 ~ Xi . 

Al ( ~vT' ~vr) Aa {-jr* ~7r ] ■ (33) 



Since the Wigner representation of number state |m) is well known [23] 

W\ m ) = (m\ Ai (xi,px) \m) 

= { —^-e-^+ti)L m [2 (xj+pl)], 
so we directly obtain the Wigner representation of L-G mode, 
w / i a ( Xl + P 2 Pi ~ ^2 ' 

' x 2 +pi p 2 - Xl 



(34) 



e-^'L mp (Q + Q 2 )L np (Qo-Q 2 ), (35) 

where Qq = p\ + p\ + x\ + x\ and Q 2 — 2p 2 x\ — 2p\x 2 . Eq. ([3"5]) is in agreement with the result of 
Ref. [5][7J. Our derivation seems economical. 
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5 The marginal distributions and fractional Fourier trans- 
form of \n, I) 

The fractional Fourier transform (FrFT) has been paid more and more attention within different 
contexts of both mathematics and physics. It is also very useful tool in Fourier optics and information 
optics. In this section, we examine the relation between the FrFT and the marginal distributions of 

W\ n ,l). 

For this purpose, we recall that the two-mode Wigner operator Ai (xi,pi) A 2 (£2,^2) = Ai (a) A 2 (/3) 
(a = (xi + ip\)/\2, P = (x\ + ipi)/v2) in entangled state representation (r| . Using the IWOP 
technique we have shown in |22j that Ai, 2 (c, 7) is just the product of two independent single-mode 
Wigner operators Ai (a) A 2 (/3) = Ai i2 (c, 7) i.e., 

A li2 (a, 7) = J k - r) (a + r\ e^'~^\ (36) 
where a = a — /?*, 7 = a + /3* and |r = t\ + ir 2 ) can be expressed in two-mode Fock space as [T2UT3"] 



exp 



~\r\ 2 +Ta\-T* 4 + 44X100), (37) 



which is the common eigenvector of X\ — X2 and Pi + P 2 , which obeys the eigenvector equations 
(X 1 - X 2 ) |t) - V2 n \t) , (Pi + P2) \t) = V2r 2 |r) . 

Performing the integration of Ai i2 (c, 7) over d 2 j (d 2 a) leads to the projection operator of the 
entangled state |t) (|£)) 

d 2 7 A li2 (a,7) = -\t) (t|| t=ct , 

7T 

|rf 2 crA 1 , 2 (cr,7) = i|^)(e|| ?=7 , (38) 

where |£) is the conjugate state of |r). Thus the marginal distributions for quantum states p in 
( T i: T 2) and (£1, £2) phase space are given by 

d 2 CrIF( ( T,7) = i(e|p|0l?=7 I 
7T 

d 2 1 W(<j, 1 ) = -(r\p\r)\ T=a , (39) 

7T 

respectively. Eq. (|39| shows that, for bipartite system, the marginal distributions can be calculated 
by evaluating the quantum average of p in (r| representations. 

Now we calculate the inner-product (r \n,l) . Note that a{a 2 , aia 2 , and J z = ^ ^ajai — c4a 2 

make up a close SU(2) Lie algebra, thus e 4 ^ Jx can be decomposed as 



e ^Jx = ^a\a 2 exp 

then we have 



- [a\ai - a 2 a 2 ) In 2 



e i<4«i ; (40) 



(r|n,Z) 



/ m p ! 



EH 



m — n — 2k 



k=0 



n„ + k 



(n p + k)\ '"^i k +i Um p -k+j)\ 



fe!(m p -A:)! ^ jl y (n p + fc-j)! 
x (r |m p - fc + j, n p + k- j) . 



(41) 



G 



Using the generating function of H m ^ n , we have 



(r' |m,n) = ^=H m , n (t'*,t') l\ 
v mini 



Substituting Ea.(|52|l into Eq.(gU leads to 



(r \n,l) = (_)«P 2 (m P -n p )/2 e -M 2 /2 

jm p \ (n p + fc)! 



! ^ 2 fc fe! (m p - fc)! 



p fc=0 
rip+fe 



Zj, ^! ( np + fc _ J)! -«mp-fc+i,n p +*-j [T , T) 



Thus the marginal distribution is 

d 2 7 VK(cr,7) 



p-kl rn I 

2™p^"p p ' 

7T Tip! 



y H) fc K + fc)! 

^ 2 k kl (m - fc)! 
fc=0 v p ' 



j! (n p + fc - j)! 



(42) 



(43) 



(44) 



Due to the presence of sum polynomial, the marginal distribution is not Gaussian. In a similar way, 
one can obtain the other marginal distribution in 7 direction. 

Before the end of this section, we mention the relation between the FrFT and marginal distribu- 
tion. In Ref. [23] we have proved that, in the context of quantum optics, the FrFT can be described 
as the matrix element of fractional operator exp[— ia(a\a\ + ajc^)] between (t'| and |/), i.e., 



F*[f{T')] = 



r) 



2 sin a 



■ exp 



i(\r'\ 



M 2 ) 



2 tan a 



2 sin a 



/ ( T ') = (t\ exp 



(4 



a 1 



where / (t') = (t'| /) . When |/) = e^ J * 



, the corresponding FrFT is 



ia. I a\ai + a\a,2 



T a [(t' \n,l)] = (r|exp 

= e -Mm P +n p ) ^ | nj ^ ; 



\m p ,n p ) 



(45) 



where we have used e _lQa i ai ai e iaa i ai = a\e la . Ea.(|45|) implies that the eigenequations of FrFT can 
also be (r \n, 1} with the eigenvalue being e - ta ( m p+ n p) ^ which is the superposition (|4Tj) of two- variable 
Hermite polynomials. In Ref. .25], we have proved that the two- variable Hermite polynomials 
(TVHP) is just the eigenfunction of the FrFT in complex form by using the IWOP technique 
and the bipartite entangled state representations. Here, we should emphasize that for any unitary 



two-mode operators U obeying the relation exp[— ia(a\a\ + a 2 a2)]U exp[ia{a\ai + a\ a 2)\ — U, the 
wave function (t| U \m p , n p ) is the eigenfunction of FrFT with the eigenvalue being e - ta ( m P+ n p) g]. 

Combing Eas. (|351) and (|3"5]) . one can obtain a simple formula connecting the FrFT and the 
marginal distribution of W\ n j\, 



J d^W (a, 7) = ~ \F a [(r — a \n, t 



(46) 



which is a— independent. Thus we can also obtain the marginal distribution by the FrFT. 
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6 L-G mode as generalized Wigner transform of Hermit- 
Gaussian modes 



In this section we shall reveal the relation between the L-G mode and the single variable Hermit- 
Gaussian(H-G) modes. Note that by taking the Fourier transformation of |r) with regard to T2 
followed by the inverse Fourier transformation, we can recover the entangled state |r) . In another 
word, It = Ti + it-i) can be decomposed into 



dxe ™V2r 2 



x + 



n 



V2/l 

x) 1 <g> 



Tl 

V2 



x- V2tiJ> , 



(47) 



in which \x) i (i = 1, 2) are the coordinate eigenvectors. Eq. (|4"T)) is called the Schmidt decomposition 
of |t) and indicates |r) is an entangled state [25] . 

Eq.(g7|) leads to 



p oc 






dxe lx 


/ —CO 








(m| 


J 




.J —OO 



V2r 



dxe 1 



~V2 



n 

V2 



V2 



|n) 



(48) 



It is interesting to notice that the integration [...] in Eq. (j48]) is similar to the single-mode Wigner 
operator, 

1 c 00 

due~ mp 



A(x,p) 



2tt 



U \ I u 

x ) ( x H 

2/ \ 2 



(49) 



and the left hand side of Ea. (l4"8"]) corresponds to TVHP mode (L-G mode). It might be expected 
that the L-G can be expressed in terms of the generalized Wigner transform (GWT). Actually, after 
making variable replacement, Eq. (|4"8"|) can be rewritten as 



(m, n \t) = 7r (m| A 



Tl T2_\ 



(-lf\n) 



If we introduce the following GWT, 

W g [fM (x,p) = (f\A(x )P )\v) 



(50) 



(51) 



which reduces to the usual Wigner transform under the condition \v) = |/), while for \v) = (—1)™ \ri) 
and I/) = |m) Eq. fBTj) becomes the right hand side of Eq. (l5ni) . which corresponds to the GWT. 

On the other hand, note that Eqs. (|2"21 and (|4"2")l . the left hand side of Eq. (|5D|) can be put into 
(without loss of the generality, letting m < n) 



(m,n \t) 



\J mini 



H m ,n {T,T*)e~ 



-|r| 2 /2 



(-l) n+m J^T* n - rn L^ m {TT*)e-^ 2 / 2 , 



).l 



(52) 



which indicates that the left hand side of Eo. tjSTfl) is just corresponding to the L-G mode, as well as 



(m I a:) 



- z H m (x) = h m (x) 



(53) 
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where H m (x) is single variable Hermite polynomial, and h m (x) just corresponds to the H-G mode, 
we have 



Thus, we can conclude that the L-G mode can be obtained by the GWT of two single-variable H-G 
modes. In addition, we should point out that the L-G mode can also be generated by windowed 
Fourier transform (which is often used in signal process) of two single- variable H-G modes by noticing 
the second line of Eq. (l47l) . 

In summary, we have endowed the Laguerre-Gaussian (LG) mode with new physical meaning in 
quantum optics, i.e., we find that it is just the wave function of the common eigenvector of the 
orbital angular momentum and the total photon number operators of 2-d oscillator in the entangled 
state representation. The common eigenvector can be obtained by using beam splitter with the 
phase difference (f> = ir/2 between the reflected and transmitted fields. With the aid of the Weyl 
ordering invariance under similar transforms, the Wigner representation of LG is directly obtained. 
It is shown that its marginal distributions can be calculated by the FrFT. In addition, L-G mode 
can also be considered as the generalized Wigner transform of Hermite Gaussian modes by using 
the Schmidt decomposition of the entangled state representation. 
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